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Abstract 



o 

Q I We review the work discussed and developed under the topic "Resummation" 

CN ■ at Working Group 2 "Multijet final states and energy flow" , of the HERA- 

bJQ. LHC Workshop. We emphasise the rôle played by HERA observables in the 

development of resummation tools via, for instance, the discovery and resum- 
mation of non-global logarithms. We describe the event-shapes subsequently 
On ■ developed for hadron colliders and présent resummed prédictions for the same 

using the automated resummation program CAESAR. We also point to on- 
^ ' going studies at HERA which can be of benefit for future measurements at 

^ ■ hadron colliders such as the LHC, specifically dijet Et and angular spectra and 

^ i the transverse momentum of the Breit current hémisphère. 
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1 Introduction 



Resummed calculations are an invaluable tool, both for the understanding of perturbative QCD dynam- 
ics at ail orders as well as for extracting, as accurately as possible, QCD parameters such as the strong 
coupling, quark masses and parton distribution functions. Thèse parameters, which cannot be directly 
computed from QCD perturbation theory itself, will be vital inputs in new physics searches at the LHC. 
(— I I Moreover, resummed expressions are also an important stepping stone to probing observable distribu- 

tions in régions where non-perturbative power corrections make a significant contribution. In this région 
one may expect a smearing of the resummed perturbative resuit with a non-perturbative function (for 
which one can adopt, for example, a renormalon-inspired model), and the resulting spectrum can be 
' confronted with data to test our understanding of non-perturbative dynamics. In ail thèse aspects, HERA 

data and observables have played an important rôle (sometimes significantly undenated in the literature) 
in furthering our knowledge, without which accurate studies of several observables at the LHC would 
simply not be possible. 

A concrète example of HERA's important rôle in this regard is the case of event shape distributions 
[1], theoretical studies of which led to the finding of non-global single-logarithmic [2] effects (discussed 
in more détail below). Prior to thèse studies it was widely believed that the HERA distributions, measured 
in the current hémisphère Breit frame, were trivially related to their e+e~ counterparts. Had such ideas, 
based on independent soft gluon émission by the hard partons, been applied directly to similar variables 
at the LHC, such as energy flows away from jets, the accuracy of theoretical prédictions would have been 
severely compromised leading almost certainly to erroneous claims and conclusions. 

Another area where HERA has played a vital rôle is in the testing of renormalon inspired models 
for power corrections, most significantly the dispersive approach [3] to 1/Q power corrections, tested 
against HERA event-shape distributions and mean-values [4]. The fact that HERA data seem to confirm 
such models , where one can think of the power corrections as arising from the émission of a gluon with 
trans verse momentum ©(Aqcd). is significant for the LHC. This is because the agreement of the renor- 
malon model with data demonstrates that the présence of initial state protons does not affect significantly 
the form of l/Q corrections. It thus sets limits on the additional non-perturbative contribution that may 
potentially be generated by the flight of struck partons through the proton cloud, which therefore does 



not appear to be significant. Once again it is accurate resummed prédictions [5] which have allowed us 
access to the non-perturbative domain hence strengthening our understanding of power corrections. 

One important aspect of resummed studies, till date, is that stringent comparisons of next-to- 
leading logarithmic resummed prédictions with data have only been carried out in cases involving ob- 
servables that vanish in the limit of two hard partons. Prominent examples reflecting the success of this 
program are provided by e+e^ 2 jet event shapes and DIS (1+1) jet event shapes as well as Drell-Yan 
vector boson trans verse momentum spectra at hadron colliders. At the LHC (and hadron colliders in gên- 
erai) one already bas two hard incoming partons and any observable dealing with final state jet production 
would take us beyond the tested two hard parton situation. Thus dijet event shapes at hadron colliders 
(discussed in détail later), which involve much more complicated considérations as far as the resum- 
mation goes, represent a situation where NLL resummations and power corrections are as yet untested. 
Bearing in mind the hadronic activity due to the underlying event at hadron colliders, it is important 
to test the picture of resummations and power corrections for thèse multiparton event shapes in cleaner 
environments. Thus LEP three-jet event shapes and similar 2 + 1 jet event shapes at HERA become 
important to study in conjunction with looking at resummation of event shapes at hadron colliders. 

Prédictions for several LEP and HERA three-jet event shapes already exist (see e.g [6] and for a 
full list of variables studied Réf. [7]) and at this workshop a prominent development presented was the 
proposai of several dijet event-shapes in hadron-hadron collisions and the resummed prédictions for their 
distributions [8]. 

Existing HERA data can also be usefuUy employed to study soft gluon radiation dynamics from 
multi-hard-parton ensembles, in the study of dijet Et and angular spectra. Thèse quantifies are somewhat 
différent from event shapes since one defines observables based on aggregate jet-momenta and angles 
rather than directly constructing them from final-state hadron momenta. Examples are the transverse en- 
ergy. Et, mismatch between the leading Et jets in dijet production and the azimuthal corrélation between 
jets once again refering to the highest Et jets in dijet production. For the former quantity there are 
no direct expérimental data as yet, but it is simply related to the dijet total rate in the région of symmetric 
Et cuts for which data does exist . For the latter quantity similarly there are direct expérimental data [9]. 
Thèse observables have smaller hadronisation corrections scaling as rather than 1/Q as for most 
event shapes. They thus offer a good opportunity to test the NLL perturbative prédictions alone without 
necessarily probing non-perturbative effects at the same time ^ . 

At this workshop developments were reported on extending existing calculations [10] for cone 
dijets, to différent jet algorithms, such as the kt algorithm, comparing to fixed order estimâtes and per- 
forming the leading order matching. Once the HERA data has been well described similar studies can 
be carried out for hadron-hadron dijets. In fact prédictions already exist for hadron-hadron dijet masses 
near threshold [11] but are not in a form conducive to direct comparisons with data containing neither 
the jet algorithms in the form actually employed in experiment, nor the matching to fixed order. How- 
ever thèse calculations provided a useful starting point for the calculations presented here, which should 
eventually lead to direct comparisons with data. 

Another area where HERA may play an important rôle is to establish whether unaccounted for 
small X effects may be significant in comparing theoretical resummations for e.g. vector boson pt spec- 
tra with expérimental data. It has been suggested that a non-perturbative intrinsic kt, growing steeply 
with X, is required to accomodate HERA data for semi-inclusive DIS processes [12]. When this obser- 
vation is extrapolated to the LHC kinematical région there is apparently significant small x broadening 
in the vector boson pt distribution. Similar effects may well arise in the case of the Higgs boson too. 
However DIS event shape studies in the Breit current hémisphère [5] apparently do not acquire such 
corrections since they are well described by conventional NLL resummations supported by dispersive 

'Although effects to do with intrinsic kt will eventually have to be accounted for similar to the case of Drell-Yan vector 
boson Pt spectra. 



power corrections [4], which are x independent ^. However there are some important caveats: 

• Unlike vector boson pt spectra, event shapes receive 1/Q hadronisation corrections unrelated to 
intrinsic kf. Thèse could mask terms originating from intrinsic kf which may yet contain the 
X dependence in question. 

• It has already been observed that including Hl data for Q < 30 GeV does spoil somewhat the 
agreement with the dispersive prédiction of universal power corrections to event shapes [5]. The 
origin of this effect could well be extra non-perturbative kt broadening related to the effects de- 
scribed above for vector boson pt. 

To get to the heart of this matter a usefui variable that has been suggested (see plenary talk by 
G. Salam at the first meeting of this workshop) is the modulus of the vector transverse momentum 
J2ieHc ^t,i of the current hémisphère in the DIS Breit frame. This quantity is simply related to the Drell- 
Yan Pt spectra and comparing theoretical prédictions, presented hère, with data from HERA should help 
to finalise whether additional small-x enhanced non-perturbative terms are needed to accomodate the 
data. We begin by first describing the results for hadron-hadron event shape variables, discussed by 
G. Salam at this workshop. Then we describe the progress in studying dijet Et and angular spectra 
(presented by M. Dasgupta and G. Corcella at the working group meetings). Finally we mention the 
results obtained thus far, for the Qt distribution of the current hémisphère and end with a look at prospects 
for continuing phenomenology at HERA, that would be of direct relevance to the LHC. 

2 Event shapes for hadron colliders 

Event shape distributions at hadron colliders, as has been the case at LEP and HERA, are important 
coUinear and infrared safe quantities, that can be used as tools for the extraction of QCD parameters, 
for instance a s, by comparing theory and data. In contrast however to more inclusive sources of the 
same information (e.g the ratio of 3 jet to 2 jet rates), event shape distributions provide a wealth of other 
information, some of which ought to be crucial in disentangling and further understanding the différent 
physics effects, relevant at hadron colliders. Thèse range from fixed-order prédictions to resummations, 
hadronisation corrections and, in conjunction with more detailed studies assesing the structure of, and 
rôle played by, the underlying event (beam fragmentation). 

Until recently there have only been limited expérimental studies of jet-shapes at hadron colliders 
[14] and no resummed theoretical prédictions for dijet shape variables at hadron colliders. Rapid récent 
developments (see Réf. [8] and références therein) in the field of perturbative resummations have now 
made theoretical estimâtes possible for a number of such distributions, introduced in [8] which we report 
on below. 

The three main theoretical developments that have led to the studies of Réf. [8] are: 

• Resummation for hadron-hadron dijet observables dépends on describing multiple soft gluon émis- 
sion from a System of four hard partons. The colour structure of the resulting soft anomalous di- 
mensions is highly non-trivial and was expUcitly computed by the Stony Brook group in a séries 
of papers (see e.g [11] and références therin). 

• The discovery of non-global observables [2]. The réalisation that standard resummation techniques 
based on angular ordering/independent-emission of soft gluons by the hard-parton ensemble, are 
not vaUd for observables that are sensitive to émissions in a limited angular range, has led to the 
introduction of observables that are made global by construction. This means that one can apply 
the technology developed by the Stony-Brook group to obtain accurate NLL prédictions for thèse 
observables, without having to resort to large A^c approximations. 

• The advent of automated resummation [15]. The development of generaUsed resummation for- 
mulae and powerful numerical methods to détermine the parameters and compute the functions 

^An exception is the jet broadening [13] but the x dependence there is of an entirely différent origin and nature. 
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Fig. 1 : Cut around the beam direction beyond rapidity r^max corresponding to the maximum rapidity reach of the detectors. 



thereof, has made it possible to study several variables at once rather than having to perform copi- 
ous, and in some cases previously unfeasible, calculations for each separate observable. 

We now discuss the différent types of variables defined and resummed in [8]. The first issue 
one has to deal with is the fact that expérimental detectors have a limited rapidity range, which can be 
modeled by a cut around the beam direction. 

This cut would then correspond to a position in rapidity of the edge of the most forward detector 
with momentum or energy resolution and the relevant values of the maximum rapidity for measurements 
is 3.5 units at the Tevatron and 5 units at the LHC. One may then worry about gluon émissions beyond 
this rapidity (i.e. inside the beam cut, see Fig.|2ll that émit softer gluons into the allowed rapidity range, 
outside the cônes depicted in Fig. |2l Such a configuration would of course render the observable non- 
global. 

To get around this potential problem, one can employ an idea suggested for 3-jet observables such 
as out-of-plane momentum flows in hadron-hadron collisions [16], which helps side-step the issue of 
non-globalness. We note that ail the observables studied here have the following functional dependence 
on a soft émission, k, coUinear to a given hard leg ^ (common to ail event shapes studied here and in 
other processes) 



where kt , i] and (f) are measured wrt a given hard leg and p represent the set of hard parton momenta 
including recoil against k while Q is the hard-scale of the process. We are particularly interested in 
émissions soft and coUinear to the beam (incoming) partons. Then an émission beyond the maximum 
detector rapidity i] > ??max corresponds to at most a contribution to the observable V ~ g^'^'^^''™'")''™'"' 
with ftinin = min(6i, 62) and bi and 62 are the values of b associated with colUnear émission near beam- 
partons 1 and 2. 

If one then choses to study the observable over a range of values such that 



then émissions more forward than ?7max do not affect the observable in the measured range of values. One 
can thus include the negligible contribution from this région and do the calculation as if the observable 
were global, ignoring the cut around the beam. Including the région beyond r/max does not alter the NLL 
resummed resuit in the suitably selected range Eq. |2l 

The price one has to pay is to limit the range of the study of the observable V, such that émissions 
beyond r/max make a negligible contribution. As we will mention later this is a more significant restriction 
for some variables compared to others (depending on the parameters a and b) but a range of study can 
always be found over which the observable can be treated as global. 

'in gênerai the values of parameters d, a, h and the function g dépend on the observable considered. For more détails and 
constraints on the varions parameters that ensure globalness and infrared and coUinear safety etc., see Réf. [15]. 




(1) 



L < (a + 6min)î? 



L = ]n l/V, 



(2) 




2.1 Global event shapes 

With the above caveat in place several variables can be safely studied (treated as global) over a wide 
range of values. An explicit example is the global transverse thrust defined as: 

= max — — , Tj_,g = 1 - T^^g , (3) 

i^i q±i 

where the thrust axis ut is defined in the plane transverse to the beam axis. The probability P{v), that 
the event shape is smaller than some value v behaves as: 



P{v) = exp 



zvr 



L = lnl/î;, (4) 



with G 12 = 2Cb + Cj, where Cb and Cj represent the total colour charges of the beam and jet 
(outgoing) partons. The above represents just the double-logarithmic contribution. The fuU resuit with 
control of up to next-to-leading single-logarithms in the exponent is considerably more complicated. It 
contains both the Stony-Brook colour évolution matrices as well as multiple émission effects (generated 
by phase-space factorisation). The automated resummation program CAESAR [15] is used to generate 
the NLL resummed resuit shown in Fig.|2l In this particular case the effect of the eut around the beam 
direction can be ignored for values T±^g > 0.15e~'^™'''=. We note that it is advisable to leave a safety 
margin between this value and the values included in measurement. 

Other global variables studied include the global thrust minor and the three jet-resolution threshold 
parameter 2/23 ■ For detailed définitions and studies of thèse variables, the reader is refered to [8]. 

We shall now proceed to look at two différent ways of defining event shapes in a given central 
région, which on its own would lead to non-globalness, and then adding terms that render them global. 



2.2 Forward suppressed observables 

Here we shall examine event shapes defined in a chosen central région C well away from the forward 
detector edges. 




Fig. 3: Figure depicting the central région marked C , containing the two iiard jets. 



First we define central _L momentum, and rapidity: 

Q±,c = '^q±i, vc = -pJ— ^« 9-L» (5) 
iec ^^fi *ec 

and an exponentially suppressed forward term, 

Then we can define an event shape in the central région which on its own would be non-global 
since we measure émissions just in C. The addition of 8^ to the event-shape renders the observable 
global as this term includes suitably the effect of émissions in the remaining région C. The exponential 
suppression of the added term reduces sensitivity to émissions in the forward région which in tum reduces 
the effect of the beam eut r?max considerably, pushing its impact to values of the observable where the 
shape cross-section is highly suppressed and thus too small to be of interest. 

The event shapes are constructed as described stepwise below: 

• Split C into two pièces: Up, Down 

• Define jet masses for each 

px,c = T^-f E 9^'' X = U,D. (7) 

Define sum and heavy-jet masses 

PSfi = PU,C + PD,C , PH,C = max{p{/,c, PD,c} • (8) 

Define global extension, with extra forward-suppressed term 

PS,£ = PS,C + Se ' PH,£ = PH,C + £c ■ (9) 

• Similarly: total and wide jet-broadenings 

Bt,£ = Bt,c + £c ' Bw,£ = Bw,c + £c ■ (10) 
At the double-log level the results assume an identical form to Eq.|4]with Gu representing a com- 
bination of total incoming (beam) and outgoing (jet) parton colour charges [8]. The full NLL resummed 
results have a substantially more complex form and results from CAESAR [15] are plotted in Fig.|4] 

'^There is considérable freedom on the choice of the central région. For instance this could be a région explicitly delimited 
in rapidity or the two hard jets themselves. 




Fig. 4: NLL resummed prédictions from CAESAR for tlie lieavy jet-mass and the wide jet-broadening witii tlie minimum jet 
transverse energy -Ex. min values of 50 and 200 GeV as shown. 



2.3 Indirectly global recoil observables 

Here we study observables that are defined exclusively in terms of partiales in the central région but 
are global. Such observables are already familiar from HERA studies. As an example, although the 
current-jet broadening wrt the photon axis of the DIS Breit frame involves only partiales that enter the 
aurrent hémisphère, the aurrent quark aaquires transverse momentum by recoil against remnant hémi- 
sphère partiales. This reaoil means that the observable is indireatly sensitive to émissions in the remnant 
hémisphère whiah makes the observables global. 

To aonstruat similar observables in the hadron-hadron ease we observe that by momentum conser- 
vation, the following relation holds : 

^q±i = -^q±i (11) 

whiah relates the sum of transverse momenta in C to that in the aomplementary région. Then the central 
partiales can be used to define a recoil term: 



Q±,c 



(12) 



which contains an indirect dependence on non-central émissions. 

Now we can define event shapes explicitly in terms of central particle momenta in C. Examples 
are the recoil jet-masses and broadenings 



Px,n = Px,c + '^±,c ) Bx,n = Bx,c + '^±,c > • • • (13) 

It is clear that since thèse observables are defined in terms of central partiales alone, the eut around 
the beam direction is not an issue here. There is however another potential problem. Due to the addition 
of the recoil term we lose direct exponentiation of the resuit in variable space. Exponentiation to NLL 
accuracy only holds in impact-parameter or b space . 

The physical effect in question here is similar to Drell-Yan Qt spectra where there are two com- 
peting mechanisms that lead to a given small Qt, Sudakov suppression of soft émissions and vectorial 
cancellation between barder émissions. Where the latter effect takes over (typically in the région where 
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Fig. 5: The recoil thmst minor as predicted by CAESAR, with a cutoff before the divergence. Only a small fraction of tlie 
cross-section is beyond tlie cutoff. 



single-logs are large asL ~ 1) we get a breakdown of the Sudakov resuit generated by CAESAR. This 
resuit is of the gênerai form: 

The resuit for recoil observables produced by CAESAR will contain a divergence in the single-log func- 
tion g2 and is eut before the divergence. Again for some variables this eut is at a position that significantly 
reduces the range of possible phenomenological studies. For other variables the divergence is at values 
of the observable that are sufficiently small so that only a few percent of the cross-section is beyond 
the cutoff. An example of the former is the recoil transverse thrust where 15% of the cross-section lies 
beyond the cut-off. For the recoil thrust minor, in contrast, the cutoff has only a moderate effect and 
much less of the cross-section is cutoff, due to the divergence in §2- 

Below we présent a table of the différent event shapes mentioned here and the impact of the two 
main limitations we discussed, the beam-cut r/max and the breakdown of resummation due to divergences 
of §2- Additionally we mention the expected impact of hadronisation corrections (not yet computed in 
fuU) on the différent observables as well as the form of the estimated contribution from the underlying 
event. 



Event-shape 


Impact of T^inax 


Resummation 
breakdown 


Underlying 
Event 


Jet hadronisation 




tolerable* 


none 


^ ^max / Q 




Tm,g 


tolerable 


none 


^ ^max / Q 


~ l/(VS7Q) 


y23 


tolerable 


none 




^ Vy23/Q* 


T±,£, PX,£ 


negligible 


none 






Bx,£ 


negligible 


none 




^ l/(V^Q) 


Tm,e 


negligible 


serious 


-i/Q 


~ 1/(VS7Q) 




negligible 


none 


-i/Q 




T±^TZ, PX,TZ 


none 


serious 


-i/Q 




T,n,-R., Bx,n 


none 


tolerable 




^ i/(V57Q) 




none 


intermediate* 







In the above the entries marked * are subject to uncertainty at présent. 



Further work is needed before the resummed expressions presented here can be compared with 



data including the matching to fixed order and computation of the power corrections for the various 
observables. This is currently in progress. 

Having discussed the hadron-hadron event shapes we now move on to describe resummed studies 
concerning dijet production at HERA which can also be straightforwardly extended to hadron-hadron 
collisons. 

3 Dijet pt and angular spectra 

It has been known for some time that dijet total rates cannot be predicted within fixed-order QCD if 
symmetric cuts are apphed to the two highest pt dijets [17]. While it was understood that the problems 
are to do with constraints on soft gluon émission, the exact nature of this constraint was only made clear 
in Réf. [10]. There it was pointed out that there are large double logarithms (aside from single logarithms 
and less singular pièces) in the slope a' (A) of the total rate, as a function of A the différence in minimum 
Pt values of the two highest pt jets. Thèse logarithms were resummed and it was shown that the slope of 
the total rate cr' as A 0. This leads to a physical behaviour of the total rate as reflected by the 
data [9]. 

To perform the comparison to data accurately however, requires two improvements to be made to 
the calculations of Réf. [10]. Firstly the exact same jet algorithm has to be employed in the theoretical 
calculations and expérimental measurements. The current algorithm used by Hl and ZEUS experiments 
is the inclusive kt algorithm. At hadron colliders variants of the cone algorithm are used and it is in 
fact a cone algorithm that was employed in Réf. [10]. However the détails of the calculation need to 
be ammended to define the cônes in ry, (/> space as is done experimentally and calculations concerning 
this were presented at the working group meeting. The second important step is matching to fixed order 
estimâtes. We report below on the leading order matching to DISENT [18] while a fuU NLO matching 
is still awaited. 

We also introduce and study two variables of related interest, the first is the différence in pt, 
between the highest pt jets ^Pt,jj = Pti — Pt2 (note that here we talk about the pt différence rather than 
the différence in the minimum Ecut, that we mentioned earlier. The resummation of this distribution 
— is essentially identical to that carried out in Réf. [10], except that here we compute the next-to- 
leading logarithms in différent versions of the jet algorithm, which should help with direct expérimental 
comparsions. We also perform the leading-order matching to DISENT. 

Having developed the calculational techniques for da/dAptjj it is then straightforward to gener- 
ate the results for the distribution in azimuthal angle between jets da/d(j)jj which requires resummation 
in the région = tt. Thèse distributions have been measured at HERA and the Tevatron (most re- 
cently by the DO collaboration). Comparing the resummation with data would represent an interesting 
challenge for the theory insofar as the status of resummation tools is concemed, and is potentially very 
instructive. 

3.1 The Aptjj and (f>jj distributions 

We shall consider dijet production in the DIS Breit frame. For the jet définition we can consider either an 
r], (p cone algorithm (such as the infrared and coUinear safe midpoint cone algorithm) or the inclusive kt 
algorithm. We shall point out to what level the two algorithms would give the same resuit and where they 
can be expected to differ. We shall use a four-vector recombination scheme where the jet four-momentum 
is the sum of individual constituent hadron four-momenta. We also impose cuts on the highest pt jets 
such that |ryi,2| < 1 and]3fi,t2 > £^mm- 

We then consider the quantity Aptjj = Pa — Pt2 which vanishes at Bom order and hence the 
distribution atthis order isjust ex: ô{ptjj). 

Beyond leading order the kinematical situation in the plane normal to the Breit axis is represented 



as before [10]: 



m = m(i,o) (15) 

Pt2 = Pt2 (cos(7r ± e), sm(7r ± e)) (16) 
kt = kt (coscpjSincp) (17) 

Thus we are considering a small déviation from the Born configuration of jets back-to-back in 
azimuth, induced by the présence of a soft gluon with trans verse momentum kt <^ Pti,t2 (which is 
not recombined by the algorithm with either hard parton) and with azimuthal angle (p. In the above e 
represents the recoil angle due to soft émission. We then have 

^Pt,jj = \Pn -Pt2\^ \ktcos(p\, (18) 

which accounts for the recoil e to first order and hence is correct to NLL accuracy. Thus for the émission 
of several soft gluons we have the pt mismatch given by 

^Pt,jj = (19) 

where k^ dénotes the single component of gluon transverse momentum, along the direction of the hard 
jets, which are nearly back-to-back in the transverse plane. The sum includes only partons not merged 
by the algorithm into the highest Et jets. 

Similarly for the dijet azimuthal angle distribution^ , we have : 

TT - 4>jj P:i —\y2kyi\. (20) 

where (pjj is the azimuthal angle between the two highest -pt jets. Note that in the above we have set pt\ = 
Pt2 = Pt since we are considering a small déviation from the Born configuration and this approximation 
is correct to NLL accuracy. We also introduced ky, the component of soft gluon momentum normal to 
the jet axis in the transverse plane. 

In either of the above two cases, i.e the Aptjj or (pjj distributions, an identical resummation is 
involved, due to the similar rôle of soft partons not recombined into jets. Henceforth we shall proceed 
with just the ^Ptjj resummation results, it being understood that similar considérations apply to (pjj in 
the région (pjj ~ vr. 

Assuming independent émission of soft gluons by the hard three-parton System (the incoming 
parton and the two outgoing partons that initiate the dijets) and factorising the phase-space Eq. as 
below^: 

(\ 1 roo f^JU 
^Pt.jj - IE^--'I = - / -j-Mb^Pt.n)Il^'"'''\ (21) 

the resummed resuit for the ^Pt,jj distribution can be expressed as 

— — — {E^^,^pt,,,)= J2 4 J2 Fa{y)C1{i,z,E^^)ws{QAPt,n). (22) 

axai^ ai^Ptjj S=q,g^ ^ ° a=T,L 

In the above ^ and z are phase-space variables that parametrise the Born dijet configuration, Fa=T,L 
dénotes the y = Q"^ /xs dependence associated to the transverse or longitudinal structure function while 
C°- is the Born matrix-element squared. The function w represents the resuit of resummation. 

^Note that the kinematical relations we dérive here would be equally valid for dijets produced in hadron-hadron collisions 
at the Tevatron or LHC and just the dynamics of multisoft gluon émission would be more complex. 

*We compute here the cross-section for the observable to be less than Aptjj from which we can easily obtain the corre- 
sponding distribution. 



The resummed expression w requires some explanation. Its form is as foUows 

roc , , 

MPt^jj) = l J Mb^PtJj) exp[-Rs{b)]S{b)qs [x/^, 1/6^ j . (23) 

Note the fact that the exponentiation holds only in b space where b is the impact parameter. The func- 
tion R{b) (we ignore the subscript S which describes either incoming quarks or gluons) is the Sudakov 
exponent which can be computed up to NLL accuracy, 

R{b) = Lg^iasL) + 52(0.^), L - ln(6Q). (24) 

while S{b) is the non-global contribution that arises from soft partons inside the jet emitting outside it. 
qs is the incoming quark or gluon density and its scale dépends on the variable b. The functions gi and 
g2 are the leading-logarithmic and next-to-leading logarithmic resummed quantities. 

For the leading logarithms gi and a subset of next-to-leading logarithms 52» generated essentially 
by exponentiation of the single-log resuit in b space, the cone and inclusive kt algorithms would give the 
same resuit, which we have computed. Starting from terms that begin with In^ b in 52 (specifically 
two soft wide-angle gluons), the foUowing two effects become important: 

• For cone algorithms the implementation of the split/merge stage affects the §2 pièce. Présent 
calculations [10] are valid to NLL accuracy if ail the energy shared by overlapping jets is given 
to the jet that would have highest pt. Note that this is différent from merging the overlapping jets 
themselves. If other merging procédures are used the calculation becomes more complex but is 
still tractable. 

• For the kt algorithm it is just being realised that running the algorithm générâtes terms that start at 
a^g \t? b in the exponent, which are not correctly treated by naive Sudakov exponentiation. Thèse 
terms, which are generated by the clustering procédure, can also be numerically accounted for in 

our case, but this is work in progress. 

The effects that we mention above cause a similar impact on the final resuit as the non-global term S{b) 
which was shown to be at around the 10% level in Réf. [10]. Hence the current results for the kt algorithm 
that do not account for the recently found additional terms and only approximately for the non-global 
logs, can be expected to change by around 10% when thèse effects will be included correctly. 

We présent in Fig. 6 preliminary results for the /^Ptjj distribution matched to the leading order 
DISENT prédiction, using the kt algorithm. The matching at présent combines quark and gluon channels 
wheras ideally one would Uke to separate the incoming quark and gluon channels with the right weights 
(0{as) coefficient functions). This would be possible if, for instance, there was parton flavour informa- 
tion explicit in the fixed order codes, a limitation of the fixed-order codes that needs to be addressed also 
for hadron-hadron event shapes to be matched to NLO prédictions. 

We also présent a comparison with HERWIG [19] results on the same quantity. The variable X in 
the figure merely refers to the effect of using the jet pt as the hard scale rather than the photon virtuahty 
Q^, formally a NNLL effect. It is amusing to note the very good agreement of the resummation with 
HERWIG but not too much can be read into it at this stage. Given the minor rôle of non-global effects 
we would expect HERWIG and our prédictions to indeed have a broad resemblence. However we should 
mention that the resummed resuit in Fig. 6 is at présent subject to change pending proper inclusion of non- 
global logs and the effect of independent soft émission at large angles. The latter is partly included in the 
results shown, through exponentiation of the one-gluon resuit as we pointed out before, but the clustering 
procédure changes this resuit at about the same level as the non-global logs {0{ag In^ b) in the exponent), 
and this feature needs to be accounted for still. Secondly the matching to LO DISENT combines channels 
and this spoils control over the In^ Q / ^Ptjj term in the expansion of the resummation to NLO. A fuU 
NLO matching with proper séparation of the channels is awaited. The HERWIG curve also includes an 
intrinsic kt component that lowers the height of the resuit at small Ptjj, which can be easily included in 
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Fig. 6: Figure showing the resummed resuit matched to fixed-order DISENT results for the variable A = QAptjj. Also 
shown, for comparison, are HERWIG results with matrix-element corrections and the DISENT resuit alone. 



the theoretical resummation but at présent is excluded. Given thèse différences the very good agreement 
one sees with HERWIG is expected to change to some extent although broadly speaking the shapes of 
the two curves are expected to be similar. Similar conclusions apply for the (pjj observable. 

4 The vector Qt of the current hémisphère 

Next we examine a quantity that, as mentioned in the introduction, makes a very good analogy with 
Drell-Yan transverse momentum, Qt, distributions. Comparison of the resummation of this observable 
with data could help to understand whether extra broadening of conventionally resummed Qt spectra, 
is generated at small x. If so this will be a significant factor at the LHC. The observable in question is 
the (modulus of) the vectorially summed transverse momenta of ail particles in the Breit frame current 
hémisphère: 

Qt = \Y.ki\- (25) 

Using momentum conservation this quantity is simply equal to the modulus of the transverse 
momenta of émissions in the remnant hémisphère. Thèse émissions can ail be ascribed to the incoming 
quark to NLL accuracy, apart from the soft wide-angle component where large-angle émissions in the 
current hémisphère can émit softer gluons into the remnant hémisphère (the by now familiar non-global 
logarithms). 

The resummed resuit for this observable can be expressed as : 

— ^ ~ ao / bdbJoibQt) exp[-R{b)]S{b)q{x, l/b^) (26) 
dQj. Jo 

where Jq is the zeroth order Bessel function, R{b) is the Sudakov exponent (the "radiator") , S{b) the 
non-global contribution and q dénotes the quark distribution summed over quark flavours with appropri- 
ate weights (charges). 

The resuit for the radiator to NLL accuracy can be expressed, as before, in terms of a leading-log 
and next-to-leading log function: 



R{b) = Lgi{a,L) + g2{a,L), L = \n{bQ). 



(27) 



We have 
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where we have A = (3oas \n[Q'^{bf], b = he"^ /2 and K = (67/18 - tï'^/Q)Ca - 5/9 n/. 

It is straightforward to express the resuit directly in Qt space and one has for the pure NLL re- 
summed terms: 
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(30) 



where R' = dR/d\n{Q/Qt 



T{l + R'/2) 

The resuit has a divergence at iî' = 2 which is due to retaining just NLL 
terms and is of the same nature as that discussed before for certain hadron-hadron event shapes and the 
Drell-Yan Qt distribution. However in the présent case the divergence is at quite low values of Qt, e.g 
for Q = 100 GeV, the divergence is at around 0.5 GeV (depending on the exact choice for Aqcd)- 
Thus it is possible to safely study the distribution down to Qt values of a few GeV using the simple 
form Eq. |3Ô1 We note that is is also possible to eliminate the divergence if one defines the radiator such 
that R{h) R{b)9(bQ — 1), which is a restriction that follows from leading-order kinematics (that one 
assumes to hold at ail orders). The résultant modification has only a negligible impact in the Qt range 
that we expect to study phenomenologically. 

After the matching to fixed-order is performed, we can probe the non-perturbative smearing e~^^^ 
that one can apply to the b space resummed resuit. Comparisons with data should hopefuUy reveal 
whether the NLL resummed resuit + 'intrinsic kt smearing, mentioned above, is sufficient at smaller 
values of x or whether extra broadening is generated in the small x région, that has a significant effect 
on the resuit. Data from Hl are already available for this distribution [20] and this should enable rapid 
developments concerning the above issue. 



5 Conclusions 

In this article we have provided a summary of the developments discussed at the HERA-LHC workshop 
working group 2, concerning the topic of all-order QCD resummations. Specifically we have mentioned 
récent work carried out for hadronic dijet event shapes, dijet Et and angular spectra and resummation of 
the current-hemisphere transverse momentum distribution in the DIS Breit frame. 

We have stressed the important rôle of HERA studies in the development of the subject from the 
LEP era and the fact that, in this regard, HERA has acted as a bridge between LEP studies of the past 
(although LEP analysis of data continues and is an important source of information) and future studies 
at both the Tevatron and the LHC. 

We have particularly tried to stress the continuing crucial rôle of HERA in testing all-order QCD 
dynamics, especially in the context of multi-hard parton observables where studies are currently ongoing. 
Careful expérimental and theoretical collaborative effort is needed here in order to confirm the picture 
developed for NLL resummations and power corrections. If this program is successful it will greatly ease 
the way for accurate QCD studies at more complex hadronic environments, such as the LHC. 
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